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ABSTRACT 

WASP-43b ("Hellie r etaUlGillon etalj) is one of the closest-orbiting hot Jupiters, with a semimajor axis a 
= 0.01526 ± 0.00018 AU and a period of only 0.81 days. However, it orbits one of the coolest stars with a 
hot Jupiter (K7V, = 4520 ± 120 K), giving the planet a modest equilibrium temperature of = 1440 ± 
40 K, assuming zero Bond albedo and uniform planetary energy redistribution. This has resulted in strong 
signal-to-noise-ratio (S/N) observations and deep eclipses in both Warm Spitzer channels (3.6 and 4.5 |J.m). 
Planets with higher S/N allow more accurate measurements of eclipse depths and brightness temperatures, 
placing tighter constraints on atmospheric composition and thermal structure. The eclipse depths and brightness 
temperatures from our jointly fit model are 0.346 ± 0.013% and 1684 ± 24 K at 3.6 ^im and 0.382 ± 0.015% 
and 1485 ± 24 K at 4.5 |J.m. The eclipse timings impr oved the estim ate of the orbital period, P, by a factor of 
three {P = 0.81347459 ± 2.1 x 10"^ days) compared to iGillon et al.l and put an upper limit on the eccentricity 
(e = 0.007^0 We use our Spitzer eclipse depths with two previously reported ground-based data points in 
the J and K bands to constrain the atmospheric properties of WASP-43b. The data rule out a strong thermal 
inversion in the dayside atmosphere of WASP-43b. Model atmospheres with no thermal inversions and fiducial 
oxygen-rich compositions are able to explain all the available data. These data, however, are insufficient to 
place stringent constraints on the molecular mixing ratios . The data suggest low day-night energy redistribution 
in the planet, consistent with previous studies, with a nominal upper limit of about 35% for the fraction of 
energy incident on the dayside, but redistributed to the nightside. 

Subject headings: Extrasolar Planets — eclipses — planets and satellites: atmospheres — planets and satellites: 
individual (WASP-43b) — techniques: photometric 



1. INTRODUCTION 



Rec ent Kepler results (Bo rucki et al.l 1201 It iBataUia et al.l 
120 121) have shown a striking increase in detections of the 
smallest candidates, and the planet-candidate lists now show 
that hot Jupiters are much less common than planets smaller 
than Neptune. But, nearly all Kepler candidates are too small, 
cold, or distant for atmospheric characterization, except the 
nearby hot Jupiters. Their host stars, bright enough for radial 
velocity (RV) measurements, subject these planets to a strong 
irradiating flux, which governs their atmospheric chemistry 
and dynamics. Thei r large sizes and large scale heights (e.g., 
IShowman & Guillot 2002) give the signal-to-noise ratio (S/N) 
needed for basic atmospheric characterization. 

The most common technique for observing hot Jupiters and 
characterizing their dayside atmospheres is secondary eclipse 
photometry (e.g., IF raine et al. 2013; Crossfield et al. 2012; 
Todorovetal. 2012; Desert et al. 2012; Deming et al.l 120 lit 
Beereretal. 2011; Demorv et al. 2007). During secondary 
eclipse, when the planet passes behind its star, we see a dip 
in integrated flux proportional to the planet-to-star flux ratio, 
or usually 0.02-0.5% in the Spitzer Space Telescope infrared 
wavelengths, where the signal is strongest. This dip is much 
lower at wavelengths accessible from the ground or from the 
Hubble Space Telescope. Techniques such as phase-curve 
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measurement (iKnutson et al.ll2009l I20T2I: iLewis et aU lMl 
ICowan et al."2012a"b'; Crossfield et al."2010'), eclipse spec- 
troscopy (Deming et al. 2013; Gibson et al. 2012; Ber ta et al.l 
2012), and ingress-egress mapping Ode Wit et al.l 120121: 
.Maieau et al..ii20 12) can reveal more than a secondary eclipse, 
but are available for only a small number of high-S/N planets. 

A secondary eclipse observed at one wavelength constrains 
an exoplanet's temperature. Multiple wavelengths constrain 
the planet's dayside spectrum, potentially yielding insight 
into the atmospheric composition and temperature structure. 
Different wavelengths probe different atmospheric levels and 
can be combined into a broadband spectr um for further 
atmospheric modelin g (e.g., |Madhusudhan & Seaged 120091; 
Stevenson e t alj|2010[) . Infrared observations are specifically 
valuable because the most abundant chemical species in plan- 
etary atmospheres (aside from H2 and He), such as H2O, 
CO, CO2, and CH4, have significant absorption and emission 
featur es at these wavelengths (e.g., iMadhusudhan & Seagerl 
'2010). Constraints on chemical composition and thermal 
structure are important for both further atmospheric modeling 
(e.g.. IShowman et aiJl2009D and studies of the planet's forma- 
tion. 

Secondary eclipse observations also provide insight into the 
exoplanet's orbit. Measuring the time of the secondary eclipse 
relative to the time of transit can establish an upper limit on 
orbital eccentricity, e, and constrain the argument of periapsis, 
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uj, independently of RV measurements. Orbital eccentricity is 
important in dynamical studies and in calculating irradiation 
levels. Apsidal precession can also be constrained by eclipse 
timing and can be used to reveal the degre e of central concen- 
tration of mass in the planetary interior (Ragozzine & Wolfl 
l2009t|Campoet al. 201 1; Lopez-Morales et al. 2010). 

WASP-43b was first detected by the Wid e- Angle Search for 
Planets (WASP) team jHellier et alj|201l !) in 2009 and 2010 
from the WASP-South and WASP-North observatories. The 
WASP team also performed follow-up measurements with the 
CORALIE spectrograph, the TRAPPIST telescope, and Eu- 
lerCAM in December 2010. These observations revealed a 
planet with a mass Mp = 1.78 Jupiter masses (Mj) and a ra- 
dius /?p = 0.93 Jupiter radii (/?j), transiting one of the coldest 
stars to host a hot Jupiter (type K7V, = 4400 ± 200 K). 
They found the planet to have an exceptionally short orbital 
period of 0.81 days and a semi-major axis of only 0.0142 AU, 
assuming the host star has a mass of = 0.58 ± 0.05 Mq. 
The planet's orbital eccentricity was constrained to e < 0.04 
at 3cr. Spectroscopic measurements of the star revealed a sur- 
face gravity of log(g) = 4.5 ± 0.2 (cgs) and a projected stellar 
rotation velocity of y, sin(0 = 4.0 ± 0.4 kms"', where / is the 
inclination of the star's pole to the line of sight. Strong Ca H 
a nd K emissio n indica tes that the star is active. 

iGillon et"an (12012 ) presented an extensive analysis of 
twenty-three transit light curves, seven occultations, and eight 
new measurements of the star's RV, observed during 2010 and 
2011 with TRAPPIST, the Very Large Telescope (VLT), and 
EulerCAM. They improved the parameters of the system sig- 
nificantly (e = 0.0035 ± 0.0043, Mp = 2.034 ± 0.052 Mj, 
= 1.036 ±0.019 Ri), refined stellar parameters (T^s = 4520 ± 
120 K, = 0.717 ± 0.025 Mq, R^ = 0.667 ± 0.01 1 Rq), and 
constrained stellar density (p* = 2.41 ± 0.08 p©). They also 
confirmed that the observed variability of the transit parame- 
ters can be attri buted to the va riability of the star itself (con- 
sistent with Hel lier et al.ll201 ll) . In addition, they detected the 
planet's thermal emission at 1.19 |.im and 2.09 |.im, inferring 
poor redistribution of heat to the night side and an atmosphere 
without a thermal inversion, using the atmospheric models of 
Ifortney et al. (2005, 2008). 

In this paper we present two secondary eclipses, observed 
with the Spitzer Space Telescope at 3.6 and 4.5 i^m using 
the Infrared Array Camera (IRAC; iFazio et al.ll2004l) . which 
further constrain the dayside emission of the planet and im- 
prove the orbital parameters of the system. We combine 
our ec lipse-depth measurements with those of ' Gillon et"al] 
(120121) to constrain the planetary atmosphere model, energy 
redistribution, and prese nce of thermal inversion using the 
iMadhusudhan & Seag ei^ retrieval method. 

The following sections present our observations; discuss 
photometric analysis; explain specific steps taken to arrive 
at the fits for each observation and a joint fit; give improved 
constraints on the orbital parameters based on available RV, 
eclipse, and transit data; discuss implications for the planetary 
emission spectrum and planetary composition; state our con- 
clusions; and, in the Appendix, supply the full set of system 
parameters from our own and previous work. The electronic 
attachment to this paper includes archival light curve files in 
FITS ASCII table and IRSA formats. 

2. OBSERVATIONS 

We observed two secondary eclipses of WASP-43b with the 
Spitzer IRAC camera in subarray mode (Program ID 70084). 
A sufficiently long basehne (Figure [U was monitored before 



the eclipses, providing good sampling of all Spitzer systemat- 
ics. To minimize intrapixel variability each target had fixed 
pointing. We used the Basic Calibrated Data (BCD) from 
Spitzer's data pipeline, version S. 18. 18.0. Basic observational 
information is given in Table [1] 



TABLE 1 
Observation Information 



Channel 


Observation 


Start Time 


Duration Exposure Number of 




Date 


[MJDuxd 


[s] Time [s] Frames 


Chi 


2011-07-30 


2455772.6845 


21421 2 10496 


Ch2 


2011-07-29 


2455771.8505 


21421 2 10496 



3. SECONDARY-ECLIPSE ANALYSIS - 
METHODOLOGY 

Exoplanet characterization requires high precision, since 
the planets' inherently weak signals are weaker than the 
systematics. In addition, Spitzer's systematics lack full 
physical characterizations. We have developed a modu- 
lar pipeline. Photometry for Orbits, Eclipses, and Transits 
(POET), that implements a wide variety of treatments of 
systematics and uses Bayesian methods to explore the pa- 
rameter space and information criteria for model choice. 
The POET pipe line is docume nted in our previous papers 
(Stevenson et al. "2010; Campo et al. 201 1'; 'Nvmever et aL 
2011; Stev enson et al . 2012; Blecic et al, 2013; Cubillos et a" 
120 13h . so we give here just a brief overview of the specific 
procedures used in this analysis. 

The pipeline uses Spitze r-supplied BCD frames to produce 
systematics-corrected light curves and parameter and uncer- 
tainty estimates, routinely achieving 85% of the photon noise 
limit or better Initially, POET masks pixels according to 
Spitzer 's permanent bad pixel masks, and then it addition- 
ally flags bad pixels (energetic particle hits, etc.) by grouping 
sets of 64 frames and performing a two-iteration, 4-a rejec- 
tion at each pixel location. Image centers with 0.01 pixel 
accuracy come from testing a variety of centering routines 
(IStevenson et al.l l20To . Supplementary Inf ormation). Sub- 
pixel 5 X interpolated aperture photometry ([Harrington et al.l 
l2007h produces the light curves. We omit frames with bad 
pixels in the photometry aperture. The background, sub- 
tracted before photometry, is an average of good pixels within 
an annulus centered on the star in each frame. 

Detector systematics vary by channel and can have both 
temporal (detector ramp) and spatial (intrapixel variability) 
components. At 3.6 and 4. 5 ixm, intrapixel sensitivi ty vari- 
ation is the dominant effect dCharbonneau et ani2005h . so ac- 
curate centering at the 0.01 pixel level is critical. We fit this 
systematic with a Bilinearly-Interpolated Su bpixel Sensitiv- 
ity (B LISS) mapping technique, following IStevenson et alJ 
(120121) . including the method to optimize the bin sizes and 
the minimum number of data points per bin. 

At 5.8, 8.0, and 16 i.im, there is temporal va riability, 
attributed to charge trapping (iKnutson et al.l l2009l) . Weak 
temporal dep enden c ies can also occur at 3.6 and 4.5 \xm 
jReach et a ll 120051: iCharbonneau et all 120051; iCampo et al] 
[2011; .Blecic et alj 120131) , while w eak spatial variability 
has been seen at 5.4 and 8.0 |xm (IStevenson et al.l [20121: 
[Anderson et al.ll2bl ll) . Thus, we consider both systematics 
in aU channels when determining the best-fit model. 
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Fig. 1. — Raw (left), binned (center, 60 points per bin, with Icr error bars), and systematics-corrected (right) secondary -eclipse light curves of WASP-43b at 
3.6 and 4.5 ^im. The results ai'e normalized to the system flux and shifted vertically for comparison. Note the different vertical scales used in each panel. The 
colored lines are best-fit models. The black curves in panel 2 ai'e models without eclipses. As seen in the binned plots of channel 2, a ramp model is not needed 
to correct for the time-dependent systematic if initial data points affected by pointing drift are clipped (see Section|3). 




mon procedure (iKnutson et al.ll20Tl b when searching for the 
best-fitting ramp. We remove the smallest number of points 
consistent with the minimal SDNR. 

Once we have found the best dataset in this way, we com- 
pare different ramp models by applying either the Akaike In- 
formation Criterion: 



10 10 10 10 

Bin Size (sec) 



10 10 10 10 

Bin Size (sec) 



AlC = x +2k, 



(2) 



Fig. 2. — Correlati ons of the residua ls for the two secondary eclipses of 
WASP-43b, following iPont et alj HWB ). The black line represents the RMS 
residual flux vs. bin size. The red line shows the predicted standard error 
scaling for Gaussian noise. The green line shows the Poisson noise. The 
black vertical lines at each bin size depict 1 cr uncertainties on the RMS resid- 
uals {RMS/ VtN, where N is the number of bins). The dotted vertical blue 
line indicates the ingress/egress timescale, and the dashed vertical green line 
indicates the eclipse duration timescale. Large excesses of several cr above 
the red line would indicate coiTelated noise at that bin size. Inclusion of 
Icr uncertainties shows no noise correlation on the timescale between the 
ingress/egress and the eclipse-duration time. 

We fit the model components simultaneously using a 
iMandel & Agoll (120021) eclipse, E(t); the time-dependent de- 
tector ramp model, R(ty, and the BLISS map, M(x,y): 



where k is the number of free parameters, or the Bayesian 
Information Criterion: 



BIC = x +'tlnA^, 



(3) 



F(x,y,t) = F,R(t)M(x,y)E(t), 



(1) 



where F(x,y,t) is the aperture photometry flux and is the 
constant system flux outside of the eclipse. 

To choose the best systematics models and aperture size, 
we analyze dozens of model comb inations and use goodness- 
of-fit criteria (Cam po et al .11201 Th . To choose the best aper- 
ture size and the number of initial points dropped during in- 
strument settling for a given channel, we minimize the stan- 
dard deviation of the normalized residuals (SDNR). Ignoring 
data points from the beginning of the observation is a com- 



where is the number of data points. The best model min- 
imizes the chosen criterion. For the final decision, the level 
of correlation in the photometric residuals is also consid- 
ered, by plotting root-mean-s quared (RMS) m o del residu- 
als vs. bin size (time interval: Po nt et af] I2006t FWinn et al.l 
l2008HCarn"po et al. 2011.) and comparing this to the theoreti- 
cal 1/v^RMS scaUng. 

We explore the phase space and estimate errors using a 
Markov-chain Monte Carlo (MCMC) routine following the 
Metropolis-Hastings random-walk algorithm, which uses in- 
dependent Gaussian proposal distributions for each parame- 
ter with widths chosen to give an acceptance rate of 30-60%. 
Each MCMC model fit begins with the Levenberg-Marquardt 
algorithm ( least-squares m inimization). We use a Bayesian 
prior (e.g., Gelmanir2002l) if a parameter's uncertainty in the 
literature is smaller than what we can a chieve. We then 
run en ough MCMC iterations to satisfy the iGelman & Rubiiil 
( 119921) convergence test. After every run, we assess conver- 
gence by examining plots of the parameter traces, pairwise 
correlations, autocorrelations, marginal posteriors, best-fitting 
model, and systematics-corrected best-fitting model. The fi- 
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nal fit is obtained from the simultaneous run of all datasets, 
sharing parameters such as the eclipse midpoint and duration 
among some or all datasets. 

We report the times of our secondary eclipses in both 
BJDjj^c (Coordinated Universal Time) and BJD^^ (BJD-j^Qg, 
Barycentric Dynamical Time), calculated using the Jet 
Propulsion Laborator y (JPL) Horizons system and following 
lEastmanet all (120101) . 

4. SECONDARY-ECLIPSE ANALYSIS - FIT DETAILS 

Light curves for both channels were extracted using every 
aperture radius from 2.00 to 4.50 pixels, in 0.25 pixel incre- 
ments. A 2D Gaussian fit found stellar centroids most con- 
sistently. We estimated the background flux using an annulus 
of 7-15 pixels from the center of the star for both channels. 
For the secondary-eclipse ingress and egress time, we used a 
Bayesian prior (f2-i = 950.5 ± 145.5 s), calculated from un- 
published WASP photometric and RV data. 

Figure [T] shows our systematics-corrected, best-fit light 
curve models. Figure|2]presents the scaling of the RMS model 
residuals vs. bin size for both channels, which shows no sig- 
nificant time correlation in the residuals. 

4. 1 . Channel 1-3.6 \im 

The most prominent systematic in this Spitzer channel is 
the intrapixel effect. The best BLISS-map bin size is 0.006 
pixels when we exclude bins with less than 4 measurements. 
The ramp and eclipse models fit without removing initial data 
points. The smallest value of BIC reveals that the best ramp 
model is quadratic; this model is 1.2 x lO"'" times more prob- 
able then the next best (linear) model. Table |2] lists the best 
ramp models, comparing their SDNR, BIC values, and eclipse 
depths. 



TABLE 2 
Channel 1 Ramp Models 



Ramp Model 


SDNR 


ABIC 


Eclipse Depth [%] 


Quadratic 


0.0039001 


0.0 


0.344 ± 0.013 


Rising 


0.0039113 


56.8 


0.292 ± 0.012 


No-Ramp 


0.0039315 


144.4 


0.268 ± 0.012 


Linear 


0.0039293 


142.1 


0.270 ± 0.012 



Figure[3]shows a comparison between the best ramp models 
and their SDNR and BIC values through all aperture sizes, 
indicating which aperture size is the best and which model 
has the lowest BIC value. 

Photometry generates consistent eclipse depths for all 
tested apertures, with the lowest SDNR at an aperture radius 
of 2.50 pixels (Figure |4|l. 

4.2. Channel 2 - 4.5 \xm 

In this channel, we noticed an upward trend in flux at the be- 
ginning of the observation, possibly due to telescope settling, 
which we do not model. We clipped 2300 initial data points 
(^38 minutes of observation), the smallest number of points 
consistent with the minimal SDNR. The 2.50 pixel aperture 
radius minimizes SDNR (Figure|5]l. 

To remove intrapixel variability, the BLISS bin size is 0.016 
pixels, ignoring bins with less than 4 points. The ramp models 
indicate a negligible ramp variation. The lowest BIC value 
corresponds to the model without a ramp (Table |3). The no- 
ramp model is 78 times more probable than the linear model. 



0.00420 

0.00413 

: 0.00406 
1 

' 0.00399 
0.00392 

0.00385- 

450 

350 

y 250 
m 

< 150 
50 
-50 




00 2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 4.50 



. No-ramp ABIC 
' Linear Ramp ABiC 
► Rising Ramp ABIC 
■ Quadr Ramp ABiC 



2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 4.50 

Photometry Aperture Size (Pixels) 

Fig. 3. — Channel 1 comparison between different ramp models. The 
plots show SDNR vs. aperture size and ABIC vs. aperture size. A lower 
SDNR value indicates a better model fit. The lowest SDNR value marks the 
best aperture size (2.50 pixels). A lower ABIC value at the best aperture 
size indicates which ramp model is the best (quadratic ramp model, green 
triangles). 
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Fig. 4. — The best-fit eclipse depths as a function of photometry aperture 
size for channel 1. The best four ramp models are plotted. The red point 
indicates the best aperture size for that channel. The eclipse-depth uncertain- 
ties are the result of 10^ MCMC iterations. The trend shows insignificant 
dependence of eclipse depth on aperture size (much less than Icr). 



0.0058 
0.0056 
Q 0.0054 
0.0052 



♦— ♦ No-ra 


Tip SDNR 




Detail 






•-• Linea 


Ramp SDNR 


0.005117 






" •-• Quad 


Ramp SDNR 


0.005116 
0.005115 










0.005114 










2 


48 2.50 


2.52 ^ — 















00 2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 4.50 



4 4 No-ramp ABIC 
" T T Linear Ramp ABIC 
. T T Quadr Ramp ABIC 




T T ▼ ^ 
T T T T 

♦ ♦ ♦ « 


T T ▼ T ▼ T ▼ 
T T ▼ T T T T 
♦ ♦♦♦♦♦♦ 



0.005a- 

50 

u 



Photometry Aperture Size (Pixels) 

Fig. 5. — Channel 2 comparison between different ramp models. Plot 
shows SDNR vs. aperture size and ABIC vs. aperture size. A lower SDNR 
value indicates a better model fit. The lowest SDNR value marks the best 
aperture size (2.50 pixels). Lower BIC values at the best aperture size indi- 
cate better models (best: no-ramp model, blue diamonds). The inset shows 
separation in SDNR for different ramp models at the best aperture size. 



We tested the dependence of eclipse depth on aperture ra- 
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TABLE 3 
Channel 2 Ramp Models 



Ramp Model 


SDNR 


ABIC 


Eclipse Depth [%] 


No-Ramp 


0.0051158 





0.392 ± 0.016 


Linear 


0.0051157 


8.8 


0.392 ±0.016 


Quadratic 


0.0051143 


14.2 


0.409 ±0.019 



dius, showing that they are all well within 1 u of each other, to 
validate the consistency of our models (Figure|6]l. 
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Fig. 6. — The best-fit eclipse depths as a function of photometry aperture 
size for channel 2. The three best ramp models are plotted. The red point 
indicates the best aperture size for that channel. The eclipse-depth uncer- 
tainties are the result of 10^ MCMC iterations. The trend shows negligible 
dependence of eclipse depth on aperture size (much less than la). 



4.3. Joint Fit 

To improve S/N on the eclipse depths, we share the eclipse 
width (duration) and eclipse midpoint phases in a joint fit of 
both datasets. The best ramp models and the best aperture 
sizes from the separate channel analyses are used in the joint 
fit. To p roduce the best joint-model fits, we iterated MCMC 
until the lGelman & Rubin (1992) diagnostics for all parame- 
ters dropped below 1%, which happened after 10^ iterations. 

The best joint-model fit parameters are in Table |4] Files 
containing the light curves, best model fits, centering data, 
photometry, etc., are included as electronic supplements to 
this article. The eclipse-midpoint time is further used for the 
subsequent orbital analysis, and the eclipse depths are used 
for the atmospheric analysis. 

5. ORBIT 

The eclipse midpoint (after a 15.2 second correction for the 
eclipse-transit light-time) has a phase of 0.5001 ± 0.0004, so 
ecoscj = 0.0001 ± 0.0006, or a 3(t upper limit of |ecoscj| < 
0.0018, consistent with a circular orbit. 

To improve the orbit solution further, we combined data 
from our observations with data from a variety of sources . 
Tra nsit midpoin t times were taken from iHellier et al.l (1201 lb 
and lGillon et alj (12012 1). and amateur observations were listed 
in the Exoplanet Transit Database (see Table ID. W e used 
CORALIE RV ob servat ions published by iHellier et al.l(l2011h 
and Gillon et al.' (120121) . No RV points needed removal due 
to the Rossiter-McLaughlin effect. We subtracted 15.2 sec- 
onds from the eclipse midpoints to correct for light-travel time 
across the orbit. We corrected all poi nts to TDB, if this was 
not already done ('Eastma n et al.l201 0'). We converted the am- 
ateur data from HJD to BJD, putting all times in a consistent 



BJDjQg format. There were 49 transit points, 23 RV points, 
and one effective eclipse obser vation. We fit all of these data 
simultaneously, as described bv lCampo et al.l(l201 II) . The free 
parameters in this fit were esinw; ecosw; P; the reference 
transit midpoint time, 7b; the RV semi-amplitude, K; and the 
RV offset, 7. The addition of the amateur transit observations 
i mproves the uncerta inty of P by a factor of three compared 
to iGillon et al.l (120121) . reducing it to 18 milliseconds. The fit 
finds an eccentricity of 0.007!j] consistent with a circular 
orbit and expectations for a close-in planet, where eccentric- 
ity should be damped by tidal interactions with the host star 
The fit results are summarized in Table |6] 

Because the bulk of the eccentricity signal for WASP-43b 
comes from esinw, and lu is practically unconstrained, it is 
likely that this eccentricity signal cor nes from the planet's in- 
teraction with the stellar tidal bulge. I Arras et aD (12012) pre- 
dict that the RV semi-amplitude of this effect is 8.9 ms 
Since our model shows that eK = ms"', it is possible 

that the majority or entirety of this signal is due to the tidal 
bulge interaction, and the true eccentricity is closer to the up- 
per limit deduced from the secondary eclipse. 

We found that, for a linear ephemeris fit, there is consid- 
erable scatter in O-C (observed time minus calculated time). 
The root-mean-square of the stated transit-time uncertainties 
is 27 seconds, while the standard deviation of the residuals is 
78 seconds. WASP-43b is close enough to its host star that 
tidal decay is a significant factor in its evolution, so we at- 
tempted to estimate the decay rate by adding a quadrati c term 
to our ephemeris model, following lAdams et al.l (120101) . Our 
model for the transit ephemeris is now: 

, A^(A^-l) 
Tm = To + PN+ 5P \ , (4) 

where is the number of orbits elapsed since the epoch To, P 
is the orbital period at Tq, and 5P = PP, where P is the short- 
term rate of change in the orbital period. Fitting this model 
to the transit data, we find that 5P = (-1.7 ± 0.3)x 10"^ days 
orbit-^ orP= -0.65 ±0.12 s year . This is illustrated in Fig- 
ure |7] This would be a 5(T detection, and is about ten times 
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Fig. 7. — O-C diagram for transit observations of WASP-43b with respect 
to the best linear ephemeris fit to the transits. Only points with an Exoplanet 
Transit Database quality rating of 3 or higher are shown here and used in 
this analysis. The red line represents the full quadratic ephemeris and shows 
that transits are progressively occuiTing earlier. The green line shows an 
ephemeris that best fits the RV and eclipse data in addition to the transit data. 
The quadratic trend depends largely on two outliers with relatively small un- 
certainties. 
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TABLE 4 

2-JoiNT Best-Fit Eclipse Light-Curve Parameters. 



Parameter 


Channel 1 


Channel 2 


Array Position (Jc, pix) 


14.99 


14.8 


Array Position (y, pix) 


15.01 


15.05 


Position Consistency" (Sx, pix) 


0.005 


0.01 


Position Consistency" ((5,., pix) 


0.006 


0.007 


Aperture Size (pix) 


2.5 


2.5 


Sky Annulus Inner Radius (pix) 


7.0 


7.0 


Sky Armulus Outer Radius (pix) 


15.0 


15.0 


System Flux Fs (pJy) 


64399.0 ±5.0 


3791 1.0 ±2.0 


Eclipse Depth (%) 


0.346 ±0.013 


0.382 ±0.015 


Brightness Temperature (K) 


1684.0 ±24.0 


1485.0 ±24.0 


Eclipse Midpoint (orbits) 


0.5002 ± 0.0004 


0.5002 ± 0.0004 


Eclipse Midpoint (BJDu^c " 2,450,000) 


5773.3172 ±0.0003 


5772.5037 ± 0.0003 


Eclipse Midpoint (BJD-j-db " 2,450,000) 


5773.3179 ±0.0003 


5772.5045 ± 0.0003 


Eclipse Duration (^4-1 , hrs) 


1.25 ±0.02 


1.25 ±0.02 


Ingress/Egress Time (^2-1 , hrs) 


0.264 ±0.018 


0.264 ±0.018 


Ramp Name 


quadramp 


no-ramp 


Ramp, Quadratic Term 


3e-05 ± 0.0004 


0± 


Ramp, Linear Term 


-0.082716 ±0.007 


0± 


Intrapixel Method 


BLISS 


BLISS 


BLISS Bin Size in x (pix) 


0.006 


0.016 


BLISS Bin Size in y (pix) 


0.006 


0.016 


Minimum Number of Points per Bin 


4 


4 


Total Frames 


10496 


10496 


Frames Used 


10124 


8004 


Rejected Frames (%) 


0.44 


1.12 


Free Parameters 


7 


2 


Number of Data Points in Fit 


10124 


8004 


AIC Value 


18137.0 


18137.0 


BIC Value 


18207.2 


18207.2 


SDNR 


0.0039006 


0.0051168 


Percentage of Photon-Limited S/N (%) 


80.3 


85.0 



'RMS frame-to-frame position difference. 



the va lue of -0.060 ± 0.015 s year"' found bv lAdams et al.l 
(I20T0I) for OGLE-TR-113b. The BIC values strongly favor 
the quadratic (decaying) ephemeris (BIC = 242) over the lin- 
ear ephemeris (BIC = 271). The transit ephemeris models are 
summarized in Table |7] 

As the circularity of WASP-43b's orbit is well established 
by the RV data and confirmed by the eclipse data for a lin- 
ear ephemeris, our secondary-eclipse point should not be ex- 
pected to deviate from an eclipse phase of 0.5 in a quadratic 
ephemeris. Assuming a circular orbit where eclipses occur at 
N+j, the eclipse time should be given by: 

T^^=To + p(^N+^-Y^-Sp(^N'-^-y (5) 

We find that our eclipse time matches this ephemeris within 
errors, lending further support to the quadratic fit. 

However, both fits look poor by eye. The standard deviation 
of the residuals about the quadratic ephemeris is 66 seconds, 
still much larger than the 27 second typical transit time uncer- 
tainty. The reduced for the quadratic fit is ^8, so much 
of the residual scatter is unexplained by both models. Ei- 
ther an actual variability in transit timing exists, or there were 
problems in the data processing or reporting. The quadratic 
trend is largely dependent on two outliers with very small un- 
certainties. Could there be mistakes in the time corrections 
for heterogeneous transit data? This is unlikely because the 
Exoplanet Transit Database indicates that all amateur obser- 
vations were submitted in UTC, while the professional data 



were unambiguous in their use of TDB. While it is possible 
that uncertainties for certain sets of transit data points may 
have been underestimated, the data come from many amateur 
and professional sources, and all would have had to make such 
errors. 

Possibly these are transit-timing variations (TTVs) due 
to interactions w ith undiscovered planets in the system 
(lAeol et al.ll2005h. an explorati on of which is beyond our 
present scope. ^Agol et al.l (12005 ') give a relation for the tran- 
sit timing deviation for a pair of planets in a mean-motion 
resonance. An Earth-sized planet in a 2:1 resonance with 
WASP-43b would produce a maximum timing deviation of 
24 seconds. Our typical deviation of 78 seconds suggests a 
lower limit of 3.2 Earth masses for such a resonance. Other 
simple-integer ratios of mean motion require larger masses to 
produce such an effect. WASP-43b's close orbit, and the tidal 
effects noted above, also raise the possibility that some of the 
signal may be due to tidal decay. A few calculations will show 
whether this is likely. 

A period change of P = -0.65 ± 0. 12 s year ' for WASP-43b 
translates to a change in semi-major axis d = (-9.1 ±1.6)x 10"^ 
AU year '.l Levrardetan(l2009l) give a relation for tidal decay, 
which for synchronous tidal rotation and negligible eccentric- 
ity and obliquity reduces to: 
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TABLE 5 
Transit Timing Data. 



Mid-Transit Time 


Uncertainty 


Source" 


Quality 


(BJDxdb) 






Rating 


2456035.66489 


0.0005 


George Hall 


2 


2456015.3273 


0.00057 


Martin Zibar 


2 


2456006.38071 


0.00101 


Frantisek Lomoz 


3 


2456006.3781 


0.00109 


Frantisek Lomoz 


3 


2456001.49859 


0.00156 


Alfonso CaiTefio 


3 


2456001.49662 


0.00035 


Gtistavo Muler Schteinman 


1 


2456001.49531 


0.00019 


Femand Entering 


1 


2455997.43508 


0.00086 


Rene Roy 


3 


2455997.43105 


0.00051 


Faustino Garcia 


3 


2455997.43008 


0.0004 


Nicolas Esseiva 


2 


2455997.42981 


0.00068 


Juanjo Gonzalez 


3 


2455984.41939 


0.00064 


Ferran Grau Hoita 


3 


2455984.41548 


0.00126 


Frantisek Lomoz 


3 


2455984.4149 


0.00047 


Fabio Martinelli 


2 


2455984.41472 


0.00071 


Frantisek Lomoz 


3 


2455979.534 


0.00044 


Nicolas Esseiva 


2 


2455979.5335 


0.0004 


Juanjo Gonzalez 


2 


2455957.57296 


0.00122 


Frantisek Lomoz 


4 


2455944.55468 


0.00106 


Roy Rene 


3 


2455940.48744 


0.0005 


Anthony Ayiomamitis 


2 


2455939.67475 


0.00052 


Ramon Naves 


2 


2455933.16473 


0.00025 


Peter Starr 


1 


2455686.68399 


0.0008 


Stan Shadick 


3 


2455682.61364 


0.00039 


Tanya Dax, Stacy Irwin 
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"The TRansiting ExoplanetS and CAndidates group (TRESCA, 
|http://var2.astro.cz/EN/tresca/index.php I supply their data to the Exo- 
planet transit Database (EID, http://var2.astro.cz/ETD/), vyhi ch per - 
forms the uniform transit analysis described by Poddanv e t al] <2010l) . 
The ETD web site provided the numbers in this table, which were con- 
verted from HJD (UTC) to BJD (TDB). 



TABLE 6 
Eccentric Orbital Model. 



Parameter 



Value 



esmuj 
ecosui'' 
e 

P (days)" 
To (BJDtdb) " 

7 (m s"')" 



-0.007 ± 0.010 
-0.0003 ± 0.0006 
007 
-87 !f 1 
P = 0.81347459 ± 2.1 xlO"' 
2,455,528.86853 ± 0.00005 
548 ±5 
-3595 ± 4 
426 



"Free parameter in MCMC fit 

where Q'^ is the ratio of the stellar tidal quality factor to 
the second-order stellar tidal Love number, k2, and is 
the stellar rotation rate. The quadratic model implies = 
2500 ±400. This is much lower than the values of 10^ - lO"' 
normally assumed. A small value of was also found by 
lAdams et alJ (l20Toh . 

Another mechanism to explain orbital decay would be the 
transfer of orbital angular momentum to another planet in the 
system instead of to the stellar rotation. If such a planet were 
in a superior orbit, then it would need a greater orbital an- 

[ ular momenturn for d ivergent planetary migration to occur 
Rodriquez et al.ll201 Ih . 

However, the inferred decay rate (regardless of cause) im- 
plies an infall time of order 10^ years, which is small com- 
pared to the estimated age of the star, so neither interpretation 



is likely. As the curves are quite divergent, fits performed af- 
ter an additional 1-2 observing seasons should resolve the 
ambiguity. 

6. ATMOSPHERE 

We modeled the dayside atmosphere of WASP-43b us- 
ing the atmospheric modehng an d retrieval method of 
iMadhusudhan & Seage J (120091 uMOi) . The model computes 
line-by-line radiative transfer in a one-dimensional, plane- 
parallel atmosphere, with constraints of local thermodynamic 
equilibrium, hydrostatic equilibrium, and global energy bal- 
ance. The pressure-temperature (P- T) profile and the molec- 
ular composition are free parameters of the model, allow- 
ing exploration of models with and without thermal inver- 
sions, and those with oxygen-rich as well as carbon-rich com- 
positions (Madhusudhan 2012). The model includes all the 
primary sources of opacity expected in hydrogen-dominated 
giant planet atmospheres in the temperature regimes of hot 
Jupiters, such as WASP-43b. The opacity sources include 
line-by-line absorption due to H2O, CO, CH4, CO2, and 
NH3, and collision-induced absorption (CIA) due to H2-H2. 
We also include hydrocarbons besides CH4, such as HCN 
and C2H2, which m ay be abundan t in carbon-rich atmo- 
spheres (Madhusudh an et alJ 1201 llA iKopparapu et alj l2012i : 
IMadhusudhan 20 1 2.) . Since in highly irradiate d oxygen-rich 
atmos pheres TiO and VO may be abundant dFortnev et alJ 
I2OO8I) . we also include line-by-line absorption due to TiO 
and VO in regions of the atmosphere where the temperatures 
exceed the corresponding condensation temperatures. Our 
molecular line data are from Freed man et al. (2008 ), Freed- 
man (personal communi c ation, 2009), Rothman et al. (2005), 
iKarkoschka & Tomask d (201(f), Karkoschka (personal com- 
munication, 2011), and Hams et al. (2008). We obtain the 
H2-H2 CIA opacities from Bor vsow et al.i(ll99 7) and B orvsowl 
(|2002). The volume mixing ratios of all the molecules are free 
parameters in the model. 

We constrain the thermal structure and composition of 
WASP-43b using our Spitzer photometric observations at 3.6 
|j.m and 4.5 |xm combined with ground-based photometric 
data using VLT/HAWK-I at 1 . 1 9 Hm and 2. 1 ^m ( Gillon et afl 
I2OI2I) . The model also has a joint constraint on the day-night 
energy redistribution and the Bond albedo, by requiring global 
energy balance - i.e., that the integrated emergent power from 
the planet does not exceed the incident irradiation. Given 
that the number of model parameters is > 10 (depending on 
the C/O ratio), and the number of available data points is 4, 
our goal is to find the regions of model space favored by the 
data, rather than to determine a unique fit. We explore the 
model parameter space using an MCMC routine (for details, 
see M adhusudhan & Seageiir2009l 120101; IMadhusudhan et alj 
1201 lai) . 

The data rule out a strong thermal inversion in the dayside 
atmosphere of WASP-43b. The data and two model spectra of 
atmospheres without thermal inversions are shown in Figure 
[8] The ground-based and Spitzer data provide complementary 
constraints on the atmospheric properties. The g round-based 
photometric bandpasses at 1.19 \im and 2.1 |a.m dGillon et alj 
I2OI2I) span spectral regions of low molecular opacity, and, 
hence, probe the deep layers of the atmosphere at pressures of 
P ^ 1 bar, beyond which the atmosphere is optically thick due 
to collision-induced opacity. Consequently, the brightness 
temperatures of the ground-based data constrain the isother- 
mal temperature structure of the deep atmosphere. On the 
other hand, the two Spitzer data sets show lower brightness 
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TABLE 7 
Ephemeris solutions 



Parameter 


Linear ephemeris (joint fit) 


Linear ephemeris (transits only) 


Quadratic ephemeris 


7b (BJDtdb) 


2455528.86853 ± 0.00005 


2455581.74441 ± 0.00004 


2455581.74433 ± 0.00004 


P (days orbit"') 


0.81347459 ± 2.1x10-'' 


0.81347395 ± 2.3 x 10"' 


0.81347788 ± 7.2x10"' 


SP (days orbit"^)" 


0.0 


0.0 


(-1.7±0.3)x 10"** 


P (s yr"')'' 


0.0 


0.0 


-0.65±0.12 






264.3 


231.7 


BIC 




271.4 


242.3 



>SP = PP, ^Derived parameter 




1 10 

A (/xm) 

Fig. 8. — Observations and model spectra for dayside thermal emission from WASP-43b. The blue fille d circles with error bars show our data in Sphzer IRAC 
channels 1 (3.6 |xm) and 2 (4.5 |j.m) and previously published ground-based data at 1.19 i^m and 2.1 |a.m jGillon et al.ll2012l) . The solid curves show the model 
spectra in the main panel, and the corresponding temperature-pressure profiles, with no thermal inversions, in the inset. The green (red) curves correspond to a 
model with solar (6 X solar) composition. Both models fit the data almost equally well. The dashed curves show blackbody spectra corresponding to planetary 
biightness temperatures of 1684 K and 1485 K, coiTesponding to the observed brightness temperatures in the Spitzer IRAC channels 1 and 2, respectively. 



temperatures at 3.6 and 4.5 relative to the ground- 
based data, which is possible only if the temperature struc- 
ture is decreasing outward in the atmosphere, causing molec- 
ular absorption in the Spitzer bands. The presence of a strong 
thermal inversion, on the contrary, would have caused molec- 
ular emission leading to higher brightness temperatures in 
the Spitzer bands relative to the ground-based bands. Con- 
sequently, the sum total of Spitzer and ground-based data rule 
out a strong thermal inversion in WASP-43b's dayside photo- 
sphere. 

The molecular composition is less well constrained by the 
data. Several physically plausible combinations of molecules 
can explain the absorption in the two Spitzer bands (e.g., 
iMadhusudhan & Seageiil2010HMadhusu dhan 2012). Figure|l] 
shows two oxygen-rich models in chemical equilibrium, with 
C/O ratio of 0.5 (solar value), but with different metallici- 
ties (solar and six times solar) and thermal profiles, both of 
which explain the data almost equally well. In both cases, 
H2O absorption in the 3.6 ^.m band, and H2O, CO, and CO2 



absorption in the 4.5 i_im band explain the Spitzer data. The 
fact that the two models fit the data equally well demonstrates 
the degeneracy between the molecular mixing ratios (via the 
metallicity) and the temperature gradient. Given the current 
photometric data, the solar metallicity model with a steep 
temperature profile (green curve) produces almost as good a 
fit as the higher metallicity model with a shallower tempera- 
ture profile (red curve), though the latter fits marginally better. 
On the other hand, a carbon-rich model with C/O >1; (e.g., 
lMadhusudhanl20 1 2b . with absorption due to CH4, CO, CjHj, 
and HCN could also explain the data. However, the current 
data are insufficient to discriminate between O-rich and C- 
rich compositions. Thus, new observations are required to 
obtain more stringent constraints on the chemical composi- 
tion of WASP-43b. Observations using the HST Wide Field 
Camera 3 in the 1.1 - 1.8 |xm bandpass can help constrain the 
H2O abundance in the atmosphere. As shown in Figure |8] an 
O-rich composition predicts strong absorption due to H2O in 
the WFC3 bandpass, which would be absent in a carbon-rich 
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atmosphere (iMadhusudhan et ani201 lal ). Similarly, observa- 
tions in other molecular bands, such as the CO band at 2.3 
\xm, can provide constraints on the corresponding molecular 
mixing ratios. 

Our observations provide nominal constraints 
on the day-night energy redistribution fraction (/r, 
IMadhusudhan & Sea ger 2009) in WASP-43b. Our model 
fits to the combined Spitzer and ground-based data allow for 
up to ^35% day-night energy redistribution in the planet. 
However, in the population of models that fit the data, models 
with higher /i values require cooler lower atmospheres on the 
dayside, and hence predict lower fluxes in the ground-based 
channels. For example, the models shown in Figure |8] have 
/r = 32 - 34%, assuming zero Bond albedo. While these 
models produce an acceptable global fit to all the four data 
points overall, they predict systematically lower fluxes in the 
ground-based channels, fitting the ground-based data points 
only at the cr lower error bars. Considering the ground- 
based points alone, without the Spitzer data, would imply a 
significantly higher continuum flux, and correspondingly a 
significantly lower day-night redistribu tion in the planet than 
--35%, consistent with the finding of lGillon et al.l(l2012h . 

The lack of a strong thermal inversion in WASP-43b is not 
surprising. At an equilibrium temperature of ^1400 K, the 
dayside atmosphere of WASP-43b is not expected to host 
gaseous TiO and V O, which have been proposed to cause 
thermal inv ersions (Spiegel et all 120091: iHubenv et al.1 [20031: 
iFortnev eta l. 2008), though hitherto unknown molecules 
that could also potentially cause such inversions cannot be 
ruled out (.Zahnle et ak 2009). The lack of a thermal inver- 
sion is also consistent with the hypothesis of iKnutson et alj 
(|2OT0l). since the host star WASP-43 is known to be active 
(lHellieretal.ll20TTI) . 



7. CONCLUSIONS 

Exoplanet secondary eclipses provide us with a unique way 
to observe the dayside spectrum of an irradiated planetary at- 
mosphere, where the opacities of the mixture of atmospheric 
trace molecules determine the thermal structure of the plane- 
tary atmosphere. 

WASP-43b has a 0.81 day period, making it one of the 
shortest-period transiting planets. It has a sm all semi-major 
axis (0.01526 ± 0.00018 AU, Gillon et al.1 120121) . The 
WA SP-43 star is a lo w mass star (M* = 0.717 ± 0.025 
Mq, iGUlon et"ani2012l) and is also one of the coldest of all 
stars hosting hot Jupiters. The close proximity of the planet 
probably induces large ti dal bulges on the planet's surface 
jRagozzine & Wolnl2009l) . The planet's projected lifetime is 
also unusually short for such a late-type host star, owing to 
tidal in-spiral. The estimated lifetime for this planet is of or- 
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der Myr (iHelheret al.ll2011l) . 

In this paper we report two Spitzer secondary eclipse ob- 
servations, using the IRAC 3.6 and 4.5 \xm channels. S/N = 
26 in channel 1 and 24 in channel 2 allowed non-ambiguous 
analysis. The final eclipse depths from our joint fit models 
are 0.346 ± 0.013% and 0.382 ± 0.015%, in channel 1 and 2 
respectively. The corresponding brightness temperatures are 
1684 ± 24 K and 1485 ± 24 K. 

Our secondary eclipse timings, along with the available RV 
data and transit photometry from the literature and amateur 
observations, provide better constraints on the orbital param- 
eters. The WASP-43b orbital period is improved by a factor of 
three (P = 0.8 1 347459 ± 2. 1 x 10"^ days). We also obtained an 
upper limit on the eccentricity. The timing of our secondary - 
eclipse observations is consistent with and suggestive of a cir- 
cular orbit with e < 0.0018. 

We combined our Spitzer eclipse depths with ground-based 
data at 1.19 and 2.1 \xm from Gillo n et al. (2012) to constrain 
the atmospheric properties of WASP-43b. The data rule out 
a strong thermal inversion in the dayside atmosphere. This 
is particularly evident from the fact that our brightness tem- 
peratures in both the Spitzer channels are lower than those ob- 
served in the ground-based channels, suggesting temperatures 
decreasing outward. Models without thermal inversions and 
with fiducial oxygen-rich compositions can explain the data. 
Current data are insufficient to provide stringent constraints 
on the chemical composition. The data do not suggest very 
efficient day-night energy redistribution in the planet, consis- 
tent with previous studies, though models with up to ^35% 
redistribution can explain the data reasonably well. 
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Thermal Emission of WASP-43b 
A. SYSTEM PARAMETERS 



TABLE 8 
System Parameters of WASP-43. 



PARAMETER 


VALUE 


REFERENCE 


Eclipse Parameters 


Eclipse Midpoint (BJDxdb) (201 1-07-30) 2455773.3179 ± 0.0003 


a 


Eclipse Midpoint (BJDxdb) (201 1-07-29) 2455772.5045 ± 0.0003 


a 


Eclipse Duration /4_| (hrs) 


1.25 ± 0.02 


a 


Eclipse Depth Spitzer IRAC, 3.6 (%) 0.346 ± 0.013 


a 


Eclipse Depth Spitzer IRAC, 4.5 (%) 0.382 ± 0.015 


a 


Eclipse Depth VLT HAWK-I, 2.095 i^m (ppm) 1560 ± 140 


c 


Eclipse Depth VLT HAWK-I, 1 . 1 86 (ppm) 790 ± 320 


c 


Ingress/Egress Time <2-i (hrs) 


0.264 ± 0.018 


a 


Orbital Parameters 


Orbital Period, P (days) 


0.81347459 ±2.1x10-' 


b 


Semimajor axis, a (AU) 


0.01526 ± 0.00018 


c 


Transit Time (BJO-ppg) 


2455726.54336 ± 0.00012 


c 


Orbital eccentricity, e 


Q Q07 +0.013 


b 


Argument of peiicenter, u) (deg) 


-87 


b 


Velocity semiamplitude. A" (ms"') 


548 ±5 


b 


Center-of-mass velocity 7 (ms'') 


-3595 ± 4 


b 


Stellar Parameters 


Spectral Type 


KV7 


c 


Mass, M* (Mq) 


0.717 ± 0.025 


c 


Radius, R* (Rq) 
Mean density, p» (P0) 


Af:f,7 +0.011 
"^■^'"-0.075 


c 
c 


Effective temperature, Teff (K) 


4520 ± 120 


c 


Surface gravity, log g* (cgs) 


4 (:4s +0.011 


c 


Projected rotation rate, v« sin(!) (km 


s'') 4.0 ± 0.4 


c 


Metallicity [Fe/H] (dex) 


-0.01 ± 0.12 


c 


Distance (pc) 


80± 20 


d 


Planetary Parameters 


Mass, Mp (Mj) 


2 034 


c 


Radius, (Rj) 


1.036 ±0.019 


c 


Surface gravity, log gp (cgs) 


3 672 +^■'"3 
j.u/i -0.012 


c 


Mean density, pp (g cm"-^) 


1 ^77 +0.063 
-0.059 


c 


Equilibiium temperature (A=0), Ta^ 


(K) 1440 


c 



a - this work (parameters derived using 2-joint fit, see Section|43) 
b - this work (see Section|5) 
c - Gillon et al. (2012) 
d - Hellier et al. (2011) 



